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Abstract 

A specific, genuinely three-dimensional mechanism of rogue wave formation, in a late stage of the 
modulational instability of a perturbed Stokes deep-water wave, is recognized through numerical 
experiments. The simulations are based on fully nonlinear equations describing weakly three- 
dimensional potential flows of an ideal fluid with a free surface in terms of conformal variables. 
Spontaneous formation of zigzag patterns for wave amplitude is observed in a nonlinear stage of 
the instability. If initial wave steepness is sufficiently high (ka > 0.06), these coherent structures 
produce rogue waves. The most tall waves appear in "turns" of the zigzags. For ka < 0.06, the 
structures decay typically without formation of steep waves. 
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The problem of rare extreme localized sea surface excitations (known as freak waves, rogue 
waves, or killer waves) has been extensively studied for last years (see, e.g., review where 
different physical mechanisms of the rogue wave phenomenon are discussed; for recent results 
in this field, see Ref . ^ , and references therein) . It has been recognized that one of the most 
important reasons for freak waves is the Benjamin- Feir (modulational) instability of a plane 
propagating wave 0, Q|. In the most probable scenario, refraction of swell in a spatially non- 
uniform current causes preliminary amplification of wave height around caustic region 
As a result, self- attractive nonlinear interaction becomes strong, so freak waves develop. 



Several simplified weakly nonlinear mode 



and its generalizations 



B, y, y, Q, 
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s, in particular the nonlinear Schroedinger equation 



12| . have been suggested to describe rogue waves 



quantitatively. But recently it has been realized that even the fully dispersive Zakharov 
equations, the most accurate among weakly nonlinear models, are still poor for describing 
rogue waves in their final stage, when the wave steepness is large 13|. It should be noted that 
a peak-to-trough height /i of a rogue wave reaches values which are typical for a corresponding 
limiting Stokes wave, that is 0.10 < h/X < 0.14, where A is a typical wave length. In the 
absence of an accurate analytical theory, there is a strong need in numerical results for 
fully nonlinear water waves. Recently, very precise simulations of exact equations for two- 
dimensional (2D) flows of an ideal fluid with a free boundary were reported and confirmed 
appearance of rogue waves in a nonlinear stage of the modulational instability Q, 1^, [ig], Q ■ 
Many of those results are based on a compact ( 1 -j- 1 )-dimensional form of governing equations 
in terms of so called conformal variables jisl . 19 . 20 . 21 . 22 1, when a shape of the boundary 
is represented parametrically a.s x = u — HY{u,t), y = Y{u,t), where x is a horizontal 
coordinate, y is a vertical coordinate, H is the Hilbert operator, and Y{u,t) is an unknown 
real- valued function. The second unknown function is a boundary value iplu, t) of the velocity 
potential. It is important that complex function Z{u, t) = u + {i — H)Y{u, t) is an analytical 
function z{w, t) evaluated at the real axis of the complex variable w = u+iv. Function z{w, t) 
corresponds to a conformal map of the lower half-plane of variable w onto the flow region 
in complex plane z = x + iy. This description effectively "straightens" the free boundary, 
and it makes possible simulations of planar potential flows on a personal computer with 
a spatial resolution up to several million points. However, real flows are always three- 
dimensional (3D), and the third spatial dimension definitely plays an important role in the 
wave dynamics. Unfortunately, the only "exact" numerical method for 3D potential flows 
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the boundary element method - is rather "expensive" despite its various sophisticated 



modifications, so the best spatial resolution availab 
points representing a free surface 
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27, 



e now is limited to a few thousand 



281 . But the problem of rogue waves 



requires much larger computational domains, at least several million points. On the other 
hand, since the modulational instability is predominantly two-dimensional, typical rogue 
waves appear long-crested. For long-crested waves, fully nonlinear approximate equations in 



30l |. The dynamical variables Y{u,q,t) 



conformal variables have been derived recently 
and q, t) now slowly vary along the second horizontal coordinate q. The corresponding 
weakly 3D theory is based on a geometric small parameter, the ratio of a typical wave 
length A to a large scale Iq along wave crests. Most recent numerical results based on 



31| . Important 3D 



this theory correspond to a spatial resolution about 16 million points 
effects have been observed. In particular, the above cited work reports observation of long- 
lived oscillating rogue waves. It has also been demonstrated there, for one example, that 
in evolution of a randomly perturbed planar wave, specific coherent zigzag structures are 
formed, which consist of obliquely oriented stripes with heightened wave amplitude. 

In the present work, 3D aspects of the Benjamin-Feir instability are considered in more 
detail. Late stages of the instability are simulated for different amplitudes of the initial 
plane wave, and then a comparison is made. The most important results of these numerical 
experiments are the following: (i) a new, genuinely 3D mechanism of freak wave formation 
has been recognized, when an extreme wave appears in a "turn" of a zigzag-shaped stripe 
of an increased wave amplitude; (ii) for the first time, two different dynamical regimes in a 
nonlinear stage of the Benjamin-Feir instability in 3D space have been observed. Let be a 
wave number, and let a be an amplitude of the first Fourier harmonics of a function Yq{u) 
representing the initial wave. If ka > 0.06 (that corresponds to ho/X > 0.02), then after 
an initial stage of the Benjamin-Feir instability, rogue waves with a ratio h/X > 0.08 are 
frequent in a wave system of a size L ~ 50A. (see FigH]). On the contrary, for ka < 0.06, 
the probability to observe a tall wave with h/X> 0.07 in a finite-size domain is very small. 
Existence of such a "threshold" is an important, essentially 3D effect. 

Equations of motion. — The nonlinear theory for surface waves on infinitely deep water, 
developed in Refs. {29! . I30I , is used herein. Equations of motion for unknown functions can be 
written in a non-canonical Hamiltonian form, involvin g va riational derivatives (5/C / Sip) and 
(SJC/SZ), where }C{iIj, Z, 2} is the kinetic energy {29, 30], and Z is the complex conjugate 



3 



of Z{u,q,t) = u + {i — H)Y{u,q,t). The Hilbert operator H is diagonal in the Fourier 
representation: it multiphes the Fourier-harmonics Ykm{t) = J Y{u,q,t)e~^^'^~'^^'^dudq by 
[i sign k] , so that 

HY{u,q,t) = j[isigJik]Ykm{t)e'^''^''^''dkdm/{2Txf. (1) 



The first equation is the so called kinematic condition on the free surface: 

\51C/5ij)' 



\Z 



u 



|2 



(2) 



The second equation is the dynamic condition (Bernoulli equation) 

'{51C/5ij)' 



ipt = -glmZ - ipuH 



\Z 



u 



|2 



Im (1 - iH) [2{6fC/6Z)Zu + {6K, / 6i))i)u] 
+ (3) 

where g is the gravitational acceleration. Equations ([2]) and ([3]) completely determine 
evolution of the system, provided the kinetic energy functional JCiip, Z, Z} is explic- 
itly given. Unfortunately, in three dimensions there is no exact compact expression for 
}C{ip, Z, Z}. However, for long-crested waves propagating mainly in x direction (the param- 
eter e ~ [X/lqY ^ 1), we have an approximate kinetic-energy functional in the form 

)C^}C = }Co + :F=-^J ipHipu du dq + JT, (4) 

where the first term /Cq describes purely 2D flows, while weak 3D corrections are introduced 
by the functional T: 

^ = ij (Zu^^, - Z,^^)G{Zu^, - Z.^u) du dq 



+ - / <! [(Z„^, - Z,^>^f/Zu] E{Z 



u 



-{Z - u) E[{Z.,^, - Z^^^Y/Z^]^dudq. (5) 

Here \1' = (1 + iH)ip, and the operators G and E are diagonal in the Fourier representation: 

Oik, m) = ^ , -, E{k, m) = —Jt}^—. (6) 

yk'^ + m? + \k\ yk'^ + m? + \k\ 

A difference between the above expression (jl]) and the unknown true water-wave kinetic 
energy functional is of order e^, since G(/c, 0) = l/{ik) for positive k, and E{k,0) = 1 (see 



Refs. [29|, |30|). Besides that, the hnear dispersion relation resulting from /C is correct in the 



entire Fourier plane. It should be noted that in Ref. [30|] a different expression for K, was 
used, with the same properties. However, if the form of /Co is fixed as given above, then the 
function G{k, m) is determined uniquely. Function E{k, m) was taken satisfying the relation 
E{k,m)/G{k,m) = ik for A; > 0. 

Numerical results and discussion. — The above equations were simulated for a randomly 
perturbed right-propagating initial wave (A = 100 m) in two sets of numerical experiments. 
The computational domain was a square 5x5 km in the first set, and 6x6 km in the second 
set, with double-periodic boundary conditions. Since initial amplitudes were taken relatively 
small {ka < 0.08, see FigHj), an exact stationary solution (Stokes wave) was approximated 
as follows: kYQ{u) ~ —0.5{kay + ka cos{ku) + (ka)"^ cos{2ku) + lA{ka)^ cos{3ku). The same 



random-phase initial perturbation was used as in Ref.[31|- The modulational instability 
developed initially in a similar way for all the runs: an exponential growth of the pertur- 
bations, and then formation of coherent zigzag structures. At the end of this first stage, 
maximal amplitudes reached values approximately 2.5 times exceeding the corresponding 
initial amplitudes (see FigH]). However, a rather sharp difference in the dynamical behav- 
ior was observed in a later stage. If ka > 0.06, then freak waves were intensively formed, 
typically in "turns" of the zigzags (see Figs. [2}|1]). Amplitudes of the freak waves exceeded 
four times the corresponding initial amplitudes, and the steepness 9 exceeded the value 7i/6 
which is limiting for a Stokes wave (not shown). On the contrary, for ka < 0.06, the waves 
remained quite smooth for most part of the time, and there was a tendency towards decay 
of the coherent structures. It is reasonable to assume that this intermediate regime should 
bring the system after some time to a weakly turbulent state, similar to what was observed 



in Refs. 



32 



331]. However, our weakly 3D model is not applicable for the weak turbulence. 



The above results imply that in a narrow parameter region ka = 0.060 ±0.005, a relative 
amplification of the initial plane wave is maximal, so freak waves are most prominent. Of 
course, the "threshold" should be understood in a statistical sense, because details of the 
dynamics depend on a particular realization of the initial random perturbation. Obviously, 
appearance of relatively high waves for ka < 0.06 cannot be excluded [see, e.g., the case 
ka = 0.055 in Fig{T](b)], but the corresponding probability is quite small, while at ka > 0.06 
the probability to have a freak wave (in a finite domain of a size about 50A) is comparable 
to 1. 
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A reason why the zigzag turns are so suitable for appearance of big waves is the following. 
It is seen in Figsj2ll3] that the wave crests in the stripes forming zigzag pattern are obliquely 
oriented with respect to x direction. Thus the group velocity has there a non-trivial q- 
component which is negative for "positively" oriented stripes (increment Ax/ Aq > along 
the stripe), and positive for "negatively" oriented stripes. Since the wave energy propagates 
basically with the group velocity, it is accumulated near the turns where increment Ax/ Aq 
changes the sign from negative to positive. This is a genuinely 3D mechanism of rogue wave 
formation, which has not been discussed earlier in the literature. Its additional numerical 
verification was performed for non-random zigzag patterns with different initial amplitudes. 
It was found that the mechanism is only actual for sufficiently high level of nonlinearity. A 
result of evolution for nearly critical initial amplitude is presented in FiglS] (compare with 
FigsHa]). 
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FIG. 1: (Color online). Maximum elevation of the free surface versus time in two sets of numerical 
experiments with a randomly perturbed planar wave of A = 27r//c = 100 m, for different initial 
amplitudes a of the first Fourier harmonics. Oscillations with a typical period about 16 s are not 
resolved in this figure. 
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FIG. 2: (Color online). Result of the Benjamin- Feir instability for ka = 0.08, t=75 min 3 s. Zigzag 
coherent structures and several big waves are seen. Here a characteristic value of the parameter e 
is approximately 0.15. 
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FIG. 3: (Color online). Sub-region of Fig 12] with a rogue wave at a zigzag turn. Wave crests in 
front of the big wave are obliquely oriented in a focusing manner. 
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FIG. 4: Profiles of the rogue wave from Figl3] 
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t=0. Surface elevation [m] 
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FIG. 5: (Color online). Formation of a freak wave from a zigzag wave envelope. Top: at t = 
maximum elevation is 1.7 m. Bottom: at t = 16 min 39.5 s maximum elevation is 5.0 m (at the 
zigzag corner). 
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